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It was claimed recently that a low order post-Newtonian (PN) Lagrangian formulation, which 
corresponds to the Euler-Lagrange equations up to an infinite PN order, can be identical to a PN 
Hamiltonian formulation at the infinite order from a theoretical point of view. This result is difficult 
to check because in most cases one does not know what both the Euler-Lagrange equations and the 
equivalent Hamiltonian are at the infinite order. However, no difficulty exists for a special IPN 
Lagrangian formulation of relativistic circular restricted three-body problem, where both the Euler- 
Lagrange equations and the equivalent Hamiltonian not only are expanded to all PN orders but 
also have converged functions. Consequently, the analytical evidence supports this claim. As far as 
numerical evidences are concerned, the Hamiltonian equivalent to the Euler-Lagrange equations for 
the lower order Lagrangian requires that they both be only at higher enough finite orders. 
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I. INTRODUCTION 

In classical mechanics, Lagrangian and Hamiltonian 
formulations are completely the same description of a 
dynamical system. Usually more attention to the Hamil¬ 
tonian formulation is paid because it has properties of a 
canonical system. 

In post-Newtonian (PN) mechanics of general relativ¬ 
ity, the two formulations are still adopted. Are they 
completely equivalent? Ten years ago two independent 
groups [1,2] answered this question. They proved the 
complete physical equivalence of the third-order post- 
Newtonian (3PN) Arnowitt-Deser-Misner (ADM) coordi¬ 
nate Hamiltonian approach to and the 3PN harmonic co¬ 
ordinate Lagrangian approach to the dynamics of spinless 
compact binaries. This result was recently extended to 
the inclusion of the next-to-next-to-leading order (4PN) 
spin-spin coupling [3]. 

However, there are two different claims on the chaotic 
behavior of compact binaries with one body spinning 
and spin effects restricted to spin-orbit (1.5PN) coupling. 
That is, the 2PN harmonic coordinate Lagrangian dy¬ 
namics allow the onset of chaos [4], but the 2PN ADM 
Hamiltonian dynamics are integrable, regular and non- 
chaotic [5,6]. 

An explanation to the opposite results was given in [7]. 
In fact, the 2PN Hamiltonian and Lagrangian formula¬ 
tions are not exactly equal but are only approximately 
related. As its detailed account, the equations of motion 
for the Lagrangian formulation use lower-order terms as 
approximations to higher-order acceleration terms in the 
Euler-Lagrange equations, while these approximations do 
not occur in the equations of motion for the Hamilto¬ 
nian formulation. It is natural that the Lagrangian has 
approximate constants of motion but the Hamiltonian 
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contains exact ones. These facts were regarded as the 
essential point for the two formulations having different 
dynamics. In this sense, the two claims that seem to be 
explicitly conflicting were thought to be correct. 

Recently, the authors of [8] revisited the equivalence 
between the Hamiltonian and Lagrangian formulations 
at PN approximations. They found that the two formu¬ 
lations at the same PN order are nonequivalent in general 
and have differences. Three simple examples of PN La¬ 
grangian formulations, including a relativistic restricted 
three-body problem with the IPN contribution from the 
circular motion of two primary objects, a spinning com¬ 
pact binary system with the Newtonian term and the 
leading-order spin-orbit coupling [8] and a binary system 
of the Newtonian term and the leading-order spin-orbit 
and spin-spin couplings [9], were used to show that the 
differences are not mainly due to the Lagrangian hav¬ 
ing the approximate Euler-Lagrange equations and the 
approximate constants of motion but come from trun¬ 
cation of higher-order PN terms between the two for¬ 
mulations transformed. An important result from the 
logic is that an equivalent Hamiltonian of a lower-order 
Lagrangian is usually at an infinite order from a theo¬ 
retical point of view or at a higher enough order from 
numerical computations. Based on this, the integrability 
or non-integrability of the Lagrangian can be known by 
that of the Hamiltonian. More recently, chaos in compa¬ 
rable mass compact binary systems with one body spin¬ 
ning was completely ruled out [10]. The reason is that a 
completely canonical higher-order Hamiltonian, which is 
equivalent to a lower-order conservative Lagrangian and 
holds four integrals of the total energy and the total an¬ 
gular momentum in an eight-dimensional phase space, 
is typically integrable [11]. This result is useful to clar¬ 
ify the doubt on the absence of chaos in the 2PN ADM 
Hamiltonian approach [5,6] and the presence of chaos in 
the 2PN harmonic coordinate Lagrangian formulation [4]. 
As a point to illustrate, two other doubts about differ¬ 
ent chaotic indicators resulting in different dynamical be¬ 
haviours of spinning compact binaries among references 
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[12-15] and different descriptions of chaotic parameter 
spaces and chaotic regions between two articles [4,16] 
have been clarified in [17-19]. 

It is worth noting that the logic result on the equiva¬ 
lence of the PN Hamiltonian and Lagrangian approaches 
at different orders is not easy to check because the exactly 
equivalent Hamiltonian of the Lagrangian is generally 
expressed as an infinite series whose convergence is un¬ 
known clearly in most cases. To provide enough evidence 
for supporting this result, we select a part of the IPN 
Lagrangian formulation of relativistic circular restricted 
three-body problem [20], where the Euler-Lagrange equa¬ 
tions can be described by a converged Taylor series and 
the equivalent Hamiltonian can also be written as an¬ 
other converged Taylor series. For our purpose, the 
Hamiltonian is derived from the Lagrangian in Sect. 2. 
Then in Sect. 3 numerical methods are used to evalu¬ 
ate whether various PN order Hamiltonians and the IPN 
Lagrangian with various PN order Euler-Lagrange equa¬ 
tions are equivalent. Finally, the main results are con¬ 
cluded in Sect. 4. 


II. POST-NEWTONIAN APPROXIMATIONS 

As in classical mechanics, a Lagrangian formulation 
L{r,r) and its Hamiltonian formulation H{p,r) satisfy 
the Legendre transformation in PN mechanics. This 
transformation is written as 

Hip,r) = p-r- L{r,r). (1) 


bodies are stationary at points Oi(a;i,0) and O 2 (x 2 , 0 ) 
with xi = —/i 2 and X 2 = /ii in the rotating reference 
frame. State variables (f, r) of the third body satisfy the 
following Lagrangian formulation 

L = Lq — kLi H—2'^2, (3) 

Lo = ^{x'^ + y'^ + x'^ + y'^) + xy - xy + U, ( 4 ) 
Li = uji{xy - xy + x'^ +y‘^), ( 5 ) 

L 2 =-^U[x‘^+y'^ + x'^+y‘^+ 2{xy - xy)]. ( 6 ) 


In the above equations, the related notations are specified 
as follows. U is of the form 



r2 


( 7 ) 


where the distances from body 3 to bodies 1 and 2 are 


ri = \J{x- xiY +y'^, 
r2 = \/{x-X 2Y +y'^. 


Lq stands for the Newtonian circular restricted three- 
body problem. Li is a IPN contribution due to the rela¬ 
tivistic effect to the circular motions of the two primaries. 
L 2 is also a IPN contribution from the relativistic effect 
to the third body, and is only a part of that in [20] for 
our purpose, wi is the IPN effect with respect to the 
angular velocity wq of the primaries and is given by 


Here r and f are coordinate and velocity, respectively. 
Canonical momentum is 


P = 


dL{r, r) 
dr 


( 2 ) 


Taking a special PN circular restricted three-body prob¬ 
lem as an example, now we derive the Hamiltonian from 
the Lagrangian in detail. 


A. Lagrangian formulation 

The circular restricted three-body problem means the 
motion of a third body (i.e. a small particle of negligible 
mass) moving around two masses mi and m 2 (mi > m 2 ). 
The two masses move in circular, coplanar orbits about 
their common center of mass, and have a constant sep¬ 
aration a and the same angular velocity. They exert a 
gravitational force on the particle but the third body 
does not affect the motion of the two massive bodies. 
Taking the unit of mass G(mi -I- m 2 ) = 1, we have the 
two masses yi = mij{mi +m 2 ) and /X 2 = m 2 /(mi -I-m 2 ). 
The unit of length requires that the constant separation 
of the two bodies should be unity. The common mean 
motion, the Newtonian angular velocity ujq, of the two 
primaries is also unity. In these unit systems, the two 


wi = (/riAi2 - 3)/(2a). (8) 

In fact, the separation a is a mark of Li and L 2 as the 
IPN effects when the velocity of light, c, is taken as one 
geometric unit in later numerical computations. 

The Lagrangian (3) is a function of velocities and co¬ 
ordinates, therefore, its equations of motion are the or¬ 
dinary Euler-Lagrange equations: 


d ,dL dL 

dt dx dx ’ dt dy dy' 

(9) 

Since the momenta Px = dL/dx and Py = dL/dy of the 
forms 

/I *^1 \ 3C/ ,. , 

Px = x-(l + —)y^ - -[x-y), 

ac^ 

(10) 

Py = 2/ + (l + — a;H- 2 ( 2 ^ + ^) 

ac^ 

(11) 

are linear functions of velocities x and y, accelerations 
can be solved exactly from Eq. (9). They have detailed 

expressions: 


Afo -I- Afi/c^ 

“ l + 317/(ac2)’ 

(12) 

To+Ti/c^ 

^ l + 3U/(ac^)' 

(13) 
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The Newtonian terms Xg and Yq and the IPN terms Xi 
and Yi are 

Xq = X + 2y + Ux, (14) 

Yo = y-2x + Uy, (15) 

Xi = 2uii{x + y) +—-L 2 -\ — [U{x + 2y) 

U a 

-{xUx + yUy){x - y)], (16) 

Yi = 2uJi{y-x) + ^L2 + ^[U{y-2x) 

-{xUx + yUy){x + y)], (17) 

where Ux = dU/dx and Uy = dU/dy. Considering that 
S = 3U/a is at the IPN level, Eqs. (12) and (13) have 
the Taylor expansions 

i ^ x„E(-i)-(|)-l + §-E(-iy(A)>|, (18) 

0 j—0 

k /c 1 

i =. y„iy;(-i)*(A)*] + h[y;(-i)i(A)i|. ( 19 ) 

i—0 j—0 

They are the Euler-Lagrange equations with PN approx- 
imations to an order A: > 1, labeled as ELk- As a point 
to illustrate, the case of A: = 0 with Xi = Yi = 0 
corresponds to the Newtonian Euler-Lagrange equations, 
marked as ELq. Erom a theoretical viewpoint, as A; —>■ 00 , 
ELk is strictly equivalent to EL given by Eqs. (12) and 
(13), namely, ELoo = EL. Note that for the generic case 
in [8], the momenta are highly nonlinear functions of ve¬ 
locities, so no exact equations of motion similar to Eqs. 
(12) and (13) but approximate equations of motion can 
be obtained from the Euler-Lagrange equations (9). This 
means that we do not know what the PN approximations 
like Eqs. (18) and (19) are converged as A; —>• 00 . 


In light of Eqs. (1), (20) and (21), we have the follow¬ 
ing Hamiltonian 

H = 2(1 ~ 

-U. (24) 

Its Taylor series at the A:th order is of the form 

k 

Hk = + 0 + ^){ypx-xpy) 

i=0 

-U. (25) 

It is clear that Hq with cui = 0 is the Newtonian Hamil¬ 
tonian formulation, and can be expressed in terms of the 
Jacobian constant Cj as Hq = —CjI2. Additionally, 
Hk is closer and closer to H as k gets larger. Without 
doubt, the exact equivalence between H and Hk should 
be Hao = H. Of course, what Hk is converged as k ^ 00 
is still unknown for the general case in [ 8 ]. 

It should be emphasized that ELk is the A;th order PN 
approximation to the Euler-Lagrange equations EL that 
is exactly derived from the IPN Lagrangian L, and Hk 
is the A;th order PN approximation to the Hamiltonian 
H. Because of the exact equivalence between EL and H, 
ELk is the A:th order PN approximation to the Hamilto¬ 
nian H, and Hk is the A;th order PN approximation to the 
Euler-Lagrange equations EL. Additionally, ELoo and 
Hoc are exactly equivalent, i.e., ELoo = EL H = Hoo- 
However, it would be up to a certain higher enough fi¬ 
nite order k rather than up to the infinite order k that 
the equivalence ELk Hk can be checked by numerical 
methods. See the following numerical investigations for 
more details. 


III. NUMERICAL INVESTIGATIONS 


B. Hamiltonian formulations 


The velocities x and y obtained from Eqs. (10) and 
( 11 ) are expressed as 


l + J/c 
Pv 

1 +J/c 


+ (1 + ^) 2 /, 


Of course, they can be expanded to the A:th order 


K ^ 

X « Px['£i-iy{-y] + il + ‘^)y, (22) 

i=0 ^ 

k r- 

y « Py[J2{-iy{^y]-{i + '^)x. ( 23 ) 


As mentioned above, Eqs. (22) and (23) are exactly iden¬ 
tical to Eqs. (20) and (21) when k 00 . 


Besides the above analytical method, a numerical 
method is used to estimate whether these PN approaches 
have constants of motion and what the accuracy of the 
constants is. Above all, we are interested in knowing 
whether these PN approaches are equivalent. 

A. Energy errors 

An eighth- and ninth-order Runge-Kutta-Fehlberg al¬ 
gorithm of variable time-steps is used to solve each of the 
above Euler-Lagrange equations ELk and Hamiltonians 
Hk. Parameters and initial conditions are Cj = 3.12, 
/i 2 = 0.001, X = y = 0.55 and i = 0. Note that the ini¬ 
tial positive value of y is given by the Jacobian constant. 
This orbit in the Newtonian problem Lq is a Kolmogorov- 
Arnold-Moser (KAM) torus on the Poincare section y = 0 
with y > 0 in Fig. 1(a), therefore, it is regular and non- 
chaotic. The integrator can give errors of the energy Hq 
for the Lagrangian Lq in the magnitude of order 10“^^ 
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or so. The long-term accumulation of energy errors is 
explicitly present in Fig. 1(b) because the integration 
scheme itself yields an artificial excitation or damping. 
If this accumulation is neglected, the energy should be 
constant. This shows that the energy iJo is actually an 
integral of the Lagrangian Lq. However, the existence 
of this excitation or damping does not make the numer¬ 
ical results unreliable during the integration time of 10^ 
due to such a high numerical accuracy. In this sense, 
not only the integrator does not necessarily use manifold 
correction methods [21-23], but also it gives true qualita¬ 
tive results as a symplectic integration algorithm [24-27] 
does. 

When the PN terms Li and L 2 are included, what 
about the accuracy of energy integrals given by the re¬ 
lated PN approximations? Let us answer this question. 
Taking the separation between the primaries, a = 31, we 
plot Fig. 2(a) in which the errors of energies of the IPN 
Euler-Lagrange equations ELi and Hamiltonian Hi are 
shown. It is worth noting that the error of energy is esti¬ 
mated by means oi A = Hi — Hi, where Hi is regarded 
as the energy of ELi at time t and Hi is the initial en¬ 
ergy. Obviously, the error for ELi is larger in about 10 
orders of magnitude than that for Hi. This result should 
be very reasonable because differences between ELi and 
Hi exist explicitly but the canonical equations are ex¬ 
actly given by the IPN Hamiltonian Hi, as shown in the 
above analytical discussions. In other words, the differ¬ 
ence between ELi and Hi is at IPN level. Of course, 
the higher the order k gets, the smaller the difference 
between ELk and Hk becomes. This is why we can see 
from Figs. 2(a) and 2(b) that the error of the 8PN Euler- 
Lagrange equations ELg, and Hamiltonian Hg is typically 
smaller than that of the IPN Euler-Lagrange equations 
ELi and Hamiltonian Hi. Without doubt, EL and H 
should be the same in the energy accuracy if no roundoff 
errors exist in Fig. 2(c). 

In addition to evaluating the accuracy of energy in¬ 
tegrals of these PN approaches, evaluating the quality 
of these PN approaches to the Euler-Lagrange equations 
EL or the Hamiltonian H is also necessary from quali¬ 
tative and quantitative numerical comparisons. See the 
following demonstrations for more information. 


B. Qualitative comparisons 


Besides the method of Poincare sections, the method 
of Lyapunov exponents is often used to detect chaos from 
order. It relates to the description of average exponential 
deviation of two nearby orbits. Based on the two-particle 
method [28], the largest Lyapunov exponent is calculated 
by 


A = lim — In 

^00 t 


d{t) 

dioy 


(26) 


where c?(0) and d{t) are distances between the two nearby 
trajectories at times 0 and t, respectively. A globally 


stable orbit is said to be regular if A = 0 but chaotic if 
A > 0. Generally speaking, it costs a long enough time to 
obtain a stabilizing value of A from the limit. Instead, a 
quicker method to find chaos is a fast Lyapunov indicator 
[29,30], defined as 


ELI = logio 


d{t) 

dioy 


(27) 


The globally stable orbit is chaotic if this indicator in¬ 
creases exponentially with time logj^Q t but ordered if this 
indicator grows polynomially. 

It can be seen clearly from the Poincare section of Fig. 
3(a) that the dynamics of EL or H in Fig. 2(c) is chaotic. 
This result is supported by the Lyapunov exponents in 
Figs. 3(b) and 3(c) and the FLIs in Fig. 3(d) and 3(e). 
What about the dynamics of these various PN approxi¬ 
mations? The key to this question can be found in Figs. 
3(b)-3(e). Here are the related details. As shown in 
Fig. 3(b), lower order PN approximations to the Euler- 
Lagrange equations EL, such as the IPN Euler-Lagrange 
equations ELi and the 4PN Euler-Lagrange equations 
ELii, are so poorer that their dynamics are regular, and 
are completely unlike the chaotic dynamics of EL. With 
increase of the PN order k, higher order PN approxima¬ 
tions to the Euler-Lagrange equations EL become better 
and better. For example, the 8PN Euler-Lagrange equa¬ 
tions ELg allows the onset of chaos, as EL does. Seen 
particularly from the evolution curve on the Lyapunov 
exponent and time, the I2PN Euler-Lagrange equations 
EL 12 seems to be very closer to EL. These results are 
also suitable for the PN Hamiltonian approximations to 
the Hamiltonian H in Fig. 3(c). When the Lyapunov 
exponents in Figs. 3(b) and 3(c) are replaced with the 
FLIs in Figs. 3(d) and 3(e), similar results can be given. 

When the separation a = 138 is instead of a = 31 in 
Fig. 3(a), an ordered KAM torus occurs. That means 
that the EL dynamics is regular and non-chaotic. In 
Figs 3(f)-3(i), lower order PN approximations such as 
ELg (or Hg) have chaotic behaviors, but higher order 
PN approximations such as EL 12 (or H 12 ) have regular 
behaviors. 

In short, the above numerical simulations seem to tell 
us that the Euler-Lagrange equations (or the Hamilto¬ 
nian approaches) at higher enough PN orders have the 
same dynamics as the Euler-Lagrange equations EL (or 
the Hamiltonian H). There is a question of whether these 
results depend on the separation a. To answer it, we fix 
the above-mentioned orbit but let a begin at 10 and end 
at 250 in increments of I. For each given value of a, the 
FLI is obtained after integration time t = 3500. In this 
way, we have dependence of FLIs on the separations a 
in several PN Lagrangian and Hamiltonian approaches, 
plotted in Fig. 4. Here 5.5 is referred as a threshold 
value of FLI for distinguishing between the regular and 
chaotic cases at this time. That is to say, an orbit is 
chaotic when its FLI is larger than threshold but ordered 
when its FLI is smaller than threshold. In light of this, 
we do not find that there are dramatic dynamical dif- 
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ferences between the Euler-Lagrange equations EL (or 
the Hamiltonian H) and the various PN approximations 
such as the IPN Hamiltonian Hi and the IPN Euler- 
Lagrange equations ELi. However, it is clearly shown in 
Table 1 that regular and chaotic domains of smaller sep¬ 
arations a in the lowest PN approaches ELi and Hi are 
explicitly different from those in EL or H. As claimed 
above, this result is of course expected. When the order 
k gets higher and higher, EL/- and Hk have smaller and 
smaller dynamical differences compared with EL or H. 
Two points are worth noting. Eirst, the same order PN 
approaches like EL 12 and H 12 (but unlike EL and H) 
are incompletely equivalent in the dynamical behaviors 
for smaller values of a. Second, all the PN approaches 
ELi, Hi, EL 12 , Hi 2 , • • •, EL and H can still have the 
same dynamics when a is larger enough. The two points 
are due to the differences among these approaches from 
the relativistic effects depending on a; smaller values of 
a result in larger relativistic effects but larger values of a 
lead to smaller relativistic effects. 


C. Quantitative comparisons 

Now we are interested in quantitative studies on the 
various PN approximations EL/- to the Hamiltonian H 
and the various PN approximations Hk to the Euler- 
Lagrange equations EL. In other words, we want to 
know how the deviation |Ar| = \Yk — vh\ between the 
position coordinate rfc for ELk (or Hk) and the posi¬ 
tion coordinate vu for H (or EL) varies with time. To 
provide some insight into the rule on the deviation with 
time, we should consider the regular dynamics in various 
PN approximations because the chaotic case gives rise to 
exponentially sensitive dependence on initial conditions. 
Eor the sake of this purpose, the parameters and initial 
conditions unlike the aforementioned ones are Cj = 2.07, 
X = 0.68 and y = 0. When a = 140 is given in Fig. 5(a), 
the curve EL is used to estimate the accuracy of numer¬ 
ical solutions between H and EL, which begins in about 
the magnitude of 10“^"^ and is in about the magnitude of 
10“^ at time t = 10000. The difference numerical solu¬ 
tions between H and ELi is rather large. With increase 
of k, ELk is soon closer to H. For instance, ELg is ba¬ 
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TABLE I: Ordered and chaotic domains of the separation a £ [10, 250] in Fig.4 


System 

Order 

Chaos 

ELi 

17, [23,25], [27,30], [34,48], 50, 52 
[58,60],[63,84],89,91,[93,99], [129,135] 
141,156,[159,173],195,196,(219,250] 

[10,16],[18,22],26,[31,33],49,51,[53,57], 
61,62,[85,88],90,92,[100,128],[136,140] 
[142,155],157,158,[174,194],[197,218] 

Hi 

16, [18,23], [26,43], [49,53], 83, 84, 194 
[88,94], [125,132], 138, 139, 152, 178 
[157,164], [166,170], 195, [219,250] 

[10,15],17,24,25,(44,48],54,55,80,81,165 
82,(85,87],[95,124],[133,137],[140,151] 
[153,156],[171,177],[179,193],[196,218] 

ELi2 

21, 22, [24,28], [32,46], 48, 137, 138 
[54,57],[60,80],85,87,(89,95], [125,131] 
[155,169], 191, 192, 193, 214, [216,250] 

[10,20], 23, 29, 30, 31, 47, [194,213] 
[49,53],58,59,[81,84],86,88,(96,124] 
[132,136], [139,154], [170,190], 215 

Hi2 

17, 21, 22, [24,28], [32,46], 48, 137, 138 
[54,57], [60,80], 85, 87, [89,95], [125,131] 
[155,169], 191, 192, 193, 214, [216,250] 

[10,16],18,19,20,23,29,30, [194,213] 
[49,53],58,59,[81,84],86,88,[96,124],31 
[132,136],[139,154],[170,190],215,47 

EL{H) [10,12], [14,16], 21, 22, [24,28], [32,46], 48 
[54,57],[60,80],85,87,[89,95],[125,131],137 
138,[155,169],191,192,193,214,[216,250] 

13, [17,20], 23, 29, 30, 31, 47, 215 
[49,53],58,59,(81,84],86,88,(96,124] 
[132,136],[139,154],[170,190],[194,213] 




FIG. 1: Panel (a) Poincare section j/ = 0 (y > 0) of an orbit with parameters Cj = 3.12 and /i 2 = 0.001 and initial conditions 
X = y — 0.55 and i = 0 in the Newtonian problem Lq. Panel (b) Energy error AE — Hq — Hq, where Ho and Ho are respectively 
energies at times t and 0. 



FIG. 2: (color online) Energy errors AE for the related PN Lagrangian formulations with the separation a = 31. Here are 
some examples to illustrate notations. For ELi, AE = Hi — Hi, where Hi is the initial energy and the energy Hi at time t 
is obtained from the solution of ELi. For Hi, AE = Hi — Hi, where the energy Hi at time t is obtained from the solution of 
Hi. For EL, AE = H — H, where H is the initial energy and the energy H at time t is obtained from the solution of EL. 
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FIG. 3: (color online) (a) Poincare section for the orbit of Fig. 1 in the PN Euler-Lagrange equations EL with the separation 
a = 31 or a = 138. Panels (b), (c), (f) and (g) relate to Lyapunov exponents A, and panels (d), (e), (h) and (i) deal with FLIs. 
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FIG. 4: (color online) Dependence of FLIs on the separation a. 
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FIG. 5: (color online) Deviation |Ar| between position solutions of the related PN Lagrangian and Hamiltonian formulations. 
Panels (a) and (c) are the deviations from H to EL, ELi {i = 1,4, 8,12). Panels (b) and (d) deal with the deviations from EL 
to H, Hi. 
















